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Abstract. We develop a theoretic description of the photogalvanic current induced by a high frequency 
radiation in asymmetric nanostructures and show that it describes well the results of numerical simulations. 
Our studies allow to understand the origin of the electronic ratchet transport in such systems and show 
that they can be used for creation of new types of detectors operating at room temperature in a terahertz 
radiation range. 
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1 Introduction 



Since the experiments of Glass et al. performed in 1974 
[T] and the theoretical explanations developed in 1976 
|2l3j it is known that an asymmetry of crystal at a mi- 
croscopic scale can lead to emergence of macroscopic sta- 
tionary directed current when the crystal is irradiated by 
an external light source. The appearance of directed cur- 
rent induced by zero-mean force of radiation in absence of 
any external static voltage has been named the photogal- 
vanic effect. This unusual phenomenon based on interplay 
of space asymmetry, relaxation and external driving had 
been mainly discussed for interaction of light with crys- 
tals [4|5j . However, this effect has a rather generic physical 
origin and the interest to it has been significantly renewed 
recently when it became clear that it may play an impor- 
tant role for transport in bio-systems where it is difficult to 
create static forces in space and where directed transport 
may be more easily generated by some oscillating param- 
eters in presence of asymmetry at a molecular level (see 
e.g. a review [B]). In this community the phenomenon be- 
came known as ratchet, following an example of pawl and 
ratchet described by Feynman and showing that a directed 
transport in asymmetric systems at a thermal equilibrium 
is forbidden by the second law of thermodynamics [7] . 

The appearance of the photogalvanic or ratchet effect 
in various systems is described in the reviews |8|9j . The 
effect has been observed with vortices in Josephson junc- 
tion arrays [10111112] . cold atoms [13], macroporous sili- 
con membranes [14j . microfluidic channels [15] and other 
systems. The great variety of systems clearly confirms a 
generic nature of the phenomenon. 



In parallel to these ratchet studies, a technological 
progress made possible to produce artificial superlattices 
in semiconductors with two-dimensional electron gas 
(2DES). The experimental studies of superlattices of an- 
tidots in a form of disks demonstrated an important con- 
tribution of periodic orbits in the transport properties of 
2DES |16ll7j . It is important to note that the interest to 
particle dynamics on a lattice of rigid disks goes back to 
the days of Galton who in far 1889 showed the appear- 
ance of statistical laws in such systems [18]. A rigorous 
mathematical description of chaos on the Galton board 
that leads to statistical laws has been given by Sinai |19j . 
The links between the chaotic dynamics, periodic orbits 
and experimental results for 2DES transport properties in 
antidot superlattices have been established in theoretical 
studies [20] . The interest to effects of microwave radiation 
on 2DES transport appeared at a relatively early stage, 
however, during a relatively long time only the case of 
antidots with a disk shape has been considered (see e.g. 
[21122) ). Due to symmetry reasons the photogalvanic effect 
is absent in such a case. 

The experimental studies of 2DES transport in asym- 
metric structures in presence of oc-driving have been 
started in [23124125] . They demonstrated the principal ex- 
istence of photogalvanic transport but no detailed analysis 
had been done for 2D structures. The quasi-lD case [H] 
has been analyzed in more detail but the ratchet transport 
in this case is rather slow due to very slow ac-driving. 
Thus, the ratchet transport in this case was linked to 
quantum tunneling effects (see [24] and discussions in [9] ) . 
On the contrary the experiments in 2D structures [23 25 
have been performed at rather high frequencies (50GHz 
and more) showing that high frequency control of ratchet 
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transport is possible in principle. Moreover, the experi- 
ment |25| demonstrated that the photogalvanic current 
exists at room temperature. Unfortunately, there was no 
further development of this interesting research line. Prob- 
ably, the absence of theoretical understanding of the phe- 
nomenon is partially responsible for this. 

The theoretical studies of ratchet transport in asym- 
metric 2DES structures in the form of semidisks Galton 
board have been started in [26] and further developed in 
|27|28j . They used extensive numerical simulations of dy- 
namical equations in combination with simple analytical 
estimates. A rigorous analytical approach, based on the 
kinetic equation, has been developed in [29] that allowed 
to solve exactly a case at a low density of asymmetric 
scatterers of a specific form. Here we combine all these 
methods that allows us to obtain a global theoretical de- 
scription of the photogalvanic effect in asymmetric nanos- 
trutures. The development of such global theory allows 
to make clear predictions for conditions under which the 
photogalvanic effect can be used for construction of room 
temperature detectors of high frequency radiation sensi- 
tive to polarization. 

The paper has the following structure: Section II gives 
a description of various models of asymmetric nanostruc- 
tures and represents simple analytical computations com- 
pared with numerical simulations; Section III describes 
analytical results based on the kinetic equation and also 
considers a general question of ratchets in dynamical sys- 
tems with or without time reversibility; the effects of mag- 
netic field on the photogalvanic current are considered in 
Section IV; discussion of the results and possible applica- 
tion of the effect are given in Section V. 

2 Model description, analytical and numerical 
results 

In our studies we consider two main types of antidots: ori- 
ented elastic semidisks and cuts (ID intervals of length 
D) which produce specular reflection from left side and 
diffusive scattering reflection from right side (see Fig. Q]). 
The cuts model has been introduced in [29] to mimic ef- 
fect of scattering on a semidisk (right diffusive side of cut 
approximately represents circular part of semidisk). The 
kinetic equation can be solved exactly in this case. For 
the semidisks Galton board we assume that the semidisks 
form a triangular lattice (see Fig.l in [28]) with R being 
a distance between disk centers and r ( i being the disk ra- 
dius. Orientation of a semidisk on (a;, y)— plane and angles 
of elastic scattering are shown in Fig. [2] For the case of 
cuts, as in [55], it is assumed that the cuts are irregularly 
distributed in space with a concentration of cuts in a unit 
area being n c (all cuts are vertical as in Fig. []}. The ki- 
netic theory J3S] also works in a case of regular lattice of 
cuts if their density is low. 

In the limit of low density of scatterers (cuts or disks) 
the scattering time r c is 

r c = l/(n c Dv) (cuts), 

t c = \/3R 2 /(8r d v) (semidisks), 1 ' 




Fig. 1. The model of oriented scatterers in a form of vertical 
cuts of length D with a concentration n c . The scattering on 
cuts is elastic from the left side and diffusive from the right 
side. The average scattering time on cuts is r c — l/(n c D Vf), 
where Vf is the Fermi velocity of 2DES. 



.'V.. 




Fig. 2. Geometry of the semidisk scatterer in (x, y)— plane 
and definition of the angle notations. A trajectory colliding 
with the semidisk at an angle <pi is scattered at an angle 0/, 
which depends on the impact parameter d; the electric field E 
is linear-polarized under angle 9 to a;— axis. 

where v is the particle velocity (for the Fermi gas v — Vp 
where Vp is the Fermi velocity of 2DES). In addition to 
scattering on antidots we assume that there is also scatter- 
ing on impurities which is characterized by the scattering 
time Tj (see [28 29 ). Between the collisions with antidots 
and impurities the electron motion is affected only by an 
electric microwave field E cos cut linear-polarized under an- 
gle 9 to x— axis (see Fig. [2]). The force acting on electron 
is F = eE and the electron velocity at time t + t s is 

2eE 

v(i + t s ) = v(t) H cos(ut + ujt s /2) sin(wr s /2), (2) 

muj 

where e,m are electron charge and mass. 

To find the velocity of stationary flow induced by mi- 
crowave oscillations we assume that without microwave 
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the equilibrium distribution of electron velocities is given 
by the Maxwell or the Fermi-Dirac distributions at tem- 
perature T. Also it is assumed that the microwave field is 
relatively weak and it only slightly perturbs the equilib- 
rium distribution. Let us start for simplicity from the case 
of the Maxwell distribution /m(v) = exp(— mw 2 /2T)/Z 
(here Z is the normalization constant). Then, according 
to @, the perturbed distribution is 



T 2 e 2 E 2 



f(v,t,T s ) = /Af(v) X 

2evE cos(tjt + ujt s /2) sin(wT s /2). 



exp[ 



Tlo 



(3) 



where we omit the velocity independent terms. 

To compute the ratchet velocity we first consider the 
cuts model. We also make certain simplifications which al- 
low to understand the physical origin of the ratchet trans- 
port keeping a more rigorous treatment for the next Sec- 
tion. Thus, we compute the averaged perturbed distribu- 
tion function. It is obtained by averaging /(v, t, r s ) over 
all initial times t and the free flight time r s . 



/c(v)=</ Pc (T)f(v,t,T s )dT s > t 



(4) 



Here p c (j~ s ) is the probability that the scattering oc- 
curs after a time r s . Its exact expression depends on the 
geometry of the sample and temperature, but for simplic- 
ity we may assume that p c {t s ) = exp(— t s /t c )/t c where r c 
is the mean scattering time. Expanding the exponent up 
to a second order in E and averaging over t, r s gives: 



/c(v) = f M (y) x 

(r c veE) 2 1 



1 



T 2 1 + (r c w) s 



(5) 



With this distribution we can compute the ratchet flow 
velocity vj >x (in x direction) considering only one scatter- 
ing on a cut. Indeed, Vf, x is the sum of two terms orig- 
inated from scatterings on the two sides of a cut Vf >x = 
v +jX —v~ iX . The contribution from the right side is v+_ x =< 
| "v | > j~K since the outgoing direction of the particle is ran- 
dom and the average depends only on the magnitude of the 
velocity before scattering |v|. On the contrary the scatter- 
ing on the left side just changes the sign of the velocity in 
the x direction: v x , as a consequence V- tX =< v x i](v x ) > 
where rj(v x ) is the step theta-function. The l/ir factor in 
the expression of v+^ x originates from the averaging over 
the random scattering angle, it is not present for v_. x 
since the scattering on the left side of the cut is elastic 
and no additional averaging is needed. If the distribution 
is isotropic these two averages are equal, keeping this fact 
in mind we may compute the contribution to v + . x , V— jX 
only from the anisotropic (non equilibrium) part of / c (v). 
This gives 



3r?e 2 E 2 



1 



r-271 



cos 



7T/2 \ 

cos 2 (0 - 9) cos(cj))d(j) 

-tt/2 J 



2V2ttto 3 T 1 + {t c uj)- 



• cos(26») 



(6) 



where 4> is the polar angle of the velocity v, 9 is the polar- 
ization angle and the coefficient comes from the average 
value of < v 4 > computed over the equilibrium distribu- 
tion f M - 

Similar computations give the ratchet velocity vf tV in 
y— direction. Here the contribution from the right cut side 
vanishes v +>y — since after a random scattering the di- 
rections with 4> and —<p have equal probability. The con- 
tribution from the left side is given by v~ tV =< v y r)(v x ) >: 



v f,v = v -,y 
\ [ cos' 

V/-7T/2 

r 2 e 2 E 2 



3r 2 e 2 E 2 



1 



2V27rm 3 r 1 + (r c w) 



\/27rTO 3 r 1 + {t c oj)' 



■ sin(20) 



(7) 



As a result Eqs. (|6l7p give the angle ip of the direc- 
tion of ratchet flow (tan^ = Vf^ y /vf tX ) as a function of 
polarization angle 9: 



tan(V>) = -2tan(20) 



(8) 



It is important to note that in the computations above 
we assumed that a scattering event definitely occurs after 
time t c . This model is rather convenient for numerical 
simulations that allows to make a comparison with the 
above theoretical estimates. Such a model corresponds to 
scattering events randomly placed in time. For a static 
random distribution of cuts in space one should compute 
the scattering probability with a transport cross section 
that gives a numerical factor 3 instead of 2 in Eq. ([5]) (see 
next Section). 

Also above it is assumed that there are no impurities. 
In their presence the result is proportional to the proba- 
bility of scattering on antidots which is equal to the ratio 
Tijiji +t c ) = t/t c due to time ergodicity. Here and below 
t = T c Ti/(r c + Ti) is the relaxation time scale determined 
by the geometrical mean of Ti and t c . Hence, in presence 
of impurities the ratchet velocity is given by 



Vf 



r 3 e 2 E 2 



1 



cos(20) \ 

2r c V2^T 1 + (ut) 2 \ 2sin(20) ) 



(9) 



where up and down terms correspond to x and y compo- 
nents of the ratchet velocity. 

The method described above can be also applied to an- 
other form of antidots, e.g. for semidisks. The main differ- 
ence from the cut model is that the direction after collision 
4>f is no longer completely determined by the impact angle 
4>i nor completely random (see Fig. [J). In fact, it is given 
by a conditional probability g{<pf\<f>i) that depends on the 
geometry of the scatterer. Due to the scatterer asymmetry, 
in general, this distribution is not invariant under y-axis 
mirror symmetry. To obtain the expression of g{<f>f\(j)i) for 
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a semidisk scatterer it is convenient to represent g via two 
parts: 



g+(<f>f\fa), fa G (-7r/2,7r/2) 
<7-(<£/l&)> & G (7r/2,37r/2) 



(10) 



Let us first compute the contribution to g from colli- 
sions with negative x impact velocity g + that correspond 
to impact angles fa E (— 7r/2, 7r/2). The scattering is elas- 
tic on all sides of the semidisk. Integrating over the impact 
parameter d (see Fig. [2]) gives the conditional probability 
g+(cj>f\fa) of scattering in the direction <pf assuming a col- 
lision with an impact angle fa. 



9+i^Afa 



2(1 



COS ( 



](</>/) 



(ii) 



Here Xw 4 — 7r,ir— i s the characteristic theta-function of 
the interval [fa — it, tt — fa]. After scattering with an im- 
pact angle fa some values of the outgoing angle 0/ are 
forbidden since the trajectory can not cross the bulk of 
the scatterer, for example fa + ir is always forbidden, the 
X function incorporates this restriction into g+(4>f\fa). 

The contribution from the trajectories with positive 
x impact velocities is given by the distribution g^{(j)t\fa) 
for impact angles in the interval fa E (tt/2,3tt/2). The 
resulting conditional probability g^(<f)f\fa) splits in two 
terms: the first comes from the collision with the straight 
edge of the semidisk and is expressed as a delta function, 
the second is related to the collisions with the curved edge 
and is similar to the expression given in Ea.lfTTj). Hence, 



-{faf | fa 



21 cos< 



1 + cos < 



-5{4> f +fa-2w) 



<t>j-4>i 



2(1 + | cos fa\) 



X[-K-0i|,|ir-<fc|](0/) 



(12) 



One can check that the distributions obtained in Eqs. pH 
T2"|) are normalized to 1. This normalization corresponds to 
the conservation of the number of particles after collision: 

ST 9+{4>f\fa)d4>f = J * g-{(f>f\fa)d(f>f = 1. 

Another effect that was not taken into account in the 
cuts model is that the scattering probability depends on 
the impact angle fa. The probability g c (fa) that the colli- 
sion occurs under the angle fa is proportional to the length 
of the segment obtained by projecting the semidisk on a 
parallel to the impact direction : 7^(1 + |cos(^j)|), this 
leads to the probability distribution 



9c 



1 + | cos fa 
2(2 + 7r) 



(13) 



The distributions given by Eqs. (|11I12I13P and the distri- 
bution function / c (v) ([5]) give the ratchet velocities via 
expressions: 



d<p f / d 2 v\v\cos(cj)f)g((i)f\fa)g c (fa)f c (v) 



Vf,i 



v f,y= / #/ / d 2 v\v\sm((j) f )g((j) f \fa)g c (fa)f c (-v) 



(14) 




Fig. 3. Dependence of the rescaled absolute value of ratchet 
velocity vf(8)/vf(6 = 0) on the polarization angle 8 for the 
cuts and semidisks models (angle is given in radians). For the 
the semidisk scatterers (black curve with circular dots) the 
lattice constant is set to R/rd = 4.5, the impurity scattering 
time is Ti ~ 15rd/vT and the Maxwell equilibrium distribu- 
tion is characterized by a temperature T. With these settings 
the numerically determined average time interval between two 
successive collisions with semidisks is r c « 0.77i~i ~ llr^/vr 
(we note that the theoretical scattering time is r c « AArd/vr 
with the thermal velocity v = vt = yf{2T/m) in Eq. (JU). In 
the simulations for the cuts model (green curve with circular 
dots) the temperature T and the time scales t c , Ti are taken 
to be the same as their numerical values in the semidisks case. 
In both cases the microwave field amplitude and frequency are 
eEr d /T = 0.4 and cjr d /V2mT = 0.316. The smooth curves 
represent the theoretical predictions of Eq. (|9]) (blue curve for 
the cuts model) and Eq. (I15|l (red line for the semidisks model). 



In the calculation of these integrals it turns out that 
the isotropic (equilibrium) term of / c (v) vanishes that cor- 
responds to the absence of the effect at equilibrium. In the 
contribution of the anisotropic term of Eq. ([5]) only the an- 
gular integrals are different from Eq. (|6l7p while the inte- 
gral on |v| is identical and leads to the same dependence 
on system parameters. As a result we obtain 



TTT 3 e 2 E 2 



1 



Vf = 



- cos(26») 

2(Hf)r c \/WTl + M 2 V sin ( 261 ) 



(15) 



To check the obtained theoretical expressions © and 
(|15p we performed numerical simulations of the cuts and 
semidisks models. The dynamical equations are solved nu- 
merically between collisions. The Maxwell equilibrium at 
temperature T is generated with the help of the Metropo- 
lis thermalization algorithm as it is described in [28] ■ The 
computation time along one trajectory is about few hun- 
dred thousands of microwave periods. The angular depen- 
dence of the ratchet velocity is shown in Figs. 13141 It is 
in a very good agreement with the obtained theoretical 
expressions both for the cuts and semidisks models. We 
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Fig. 4. Dependence of the angle tp of the ratchet flow direction 
on the polarization angle 8 for the parameters of Fig. [2] with 
the same choice of colors: black curve for numerical data and 
red line for theory (t/j = n — 29) in the semidisks model; green 
curve for numerical data and blue curve for theory (see Eq. ((8|) 
in the cuts model; angles are given in radians. 



will consider the dependence on the parameters r c ,Ti in 
the next Section. 



3 The kinetic equation approach 

The approach described above can be also used for the 
Fcrmi-Dirac equilibrium distribution. However, it if more 
convenient to use a more general approach based on the 
kinetic equation which reads: 



df df f-fo 

— + eEcosujt^- = — 

at op t 



(16) 



where fo is an unperturbed equilibrium distribution, p = 
rav. The solution can be presented in a form of expansion 
over powers of external weak field: / = /o + /i + /a + 
The first term is 

MP, t) = - 9n T l Ev S IT exp ^) + CC ' 
2(1 + lujT) Oe 

where CC is a complex conjugated part. The time aver- 
aged correction of fi gives 

< fzfat) > t = -reE < cos Lot^p- > 

Op 



(reEv) 2 d 2 f 
2(l + (wr) 2 ) de 2 ' 



(18) 



where we use an approximation that the isotropic term 
originating from the term dfx/dv does not contribute to 
the ratchet velocity and therefore can be omitted. Also 
we assume that the relaxation time r is independent of 



particle energy. The correction (fT8"|) has the same form as 
in Eq. (O. Thus, for example, 



«/,< 



< \v\ > /ir- < v x r/(v x ) >- 



(emrE) 2 
6(l + (wr) 2 ) 



cos(26>) 



4d 2 f (v) 



de 2 



dv 



and for the Fermi-Dirac distribution at T <C Ep we obtain 
the ratchet flow velocity for the cuts model 



Vf 



T^e'E 2 1 - ir 2 T 2 /{2AE 2 F ) f - cos(29) 

2sin(26>) 



(19) 



2TTT c y / 2m 3 E F 1 + H 2 
while for the semidisks model 

^e 2 E 2 1 - 7 r 2 r 2 /(24ii; 2 ? ) /-cos(2#)\ 
2(2 + ir)T c ^2m i E F 1 + (wr) 2 ^ sin(26») ) 



Vf 



(20) 



The angle dependence remains the same as for the Maxwell 
equilibrium. 

The case of weak asymmetry can be also treated in a 
more formal way based of the Green function formalism in 
the kinetic equation. First we note that in the second order 
on alternating electric field the components of steady-state 
current density can be described by phenomenological ex- 
pressions 

jx — ^xxx\E x \ ~t~ Q?xyy\Ey\ , 

j v = Re(a yxy )(E x E* y + E*E y ) + ImK^EE*], .(21) 

The components of the photogalvanic tensor a xxx , a xyy 
and He(a yxy ) determine the response to the linear-polarized 
microwave field. For linear polarization along x or y axes 
the current flows along the x direction; the current in y 
direction appears for a tilted linear-polarized electric field. 
The kinetic equation written in the operator form reads 



df 
dt 



Ff - //, 



(22) 



where f(p, <j)) is the distribution function and p = 
p(cos(f>, sin0) is the electron momentum. The term (Ff) 
represents the action of electric field E(t) = Re(E w e~ tcJt ) 
of the electromagnetic wave with the complex amplitude 
E„ = EI W : 

„. . d , . d sin0 d . 

F = eE i - = e [E x cos0- ) + 

op op p oq> 

+Ev (^ < l > ± + ^±)] = lF u e^+CC (23) 
op p dtp 2 

The term If is the collision integral 



if( P , <t>) = / dcf>'[w(cf>\ <j>)f(p, 0') - w{<f>, <t>')f{p, <!>)], 

Jo 

(24) 

where W(4>',cj>) is the scattering probability. 
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The stationary solution of the kinetic equation in the 
second order of electric field can be written as 



f 2 = ±Re(&>F u &'F- u f 



(25) 



where F^ = eE w ^, G is the Green function of kinetic 
equation. For small asymmetry 



G" = 



1 



1 



iu>—1/t (iuj—I/t)'- 



(26) 



where I~ is the antisymmetric part of scattering operator 
on cuts, semidisks etc. The quantity r = TiT c j (r, + t c ) can 
be attributed to the symmetric part of scattering on cuts 
and to the impurity relaxation time Tj. 
The photogalvanic current reads 



iJrc / M^ E -4 } i^ (E " v)/ ° 



-r(E_ u 



d 



dp (1 — iur) 2 



(27) 



Here v = p/m is the electron velocity, fo(e) is the equi- 
librium distribution function in energy e, prime means 
derivative over energy. 

From Eq. (|27|) we find the current induced by the linear 
polarized microwave field 



ji = j-2~Re J dp^T 2 vJ VjVki^Y 



(1 — iujt)'- 



;ViVjI v k \E~E k . 



(28) 



We use the identities I£(e) = 0, < Ix(p) >= 0, and 
< pipjlpk >= — < PklpiPj >, where £(e) is an arbitrary 
function of energy, x(p) is an arbitrary function of mo- 
mentum; angular brackets stand for the operation of av- 
erage over angles in p-space: < ... >= J ^ (...). The first 
and the second identities follow from the conservation of 
number of particles, the third is the consequence of the de- 
tailed balance principle according to which the probability 
of transition has the symmetry W(p', p) = W(— p, — p'). 

Further we consider an algebraic energy dependencies 
of relaxations times: t, t c cx e s . This dependence corre- 
sponds to the scattering on a geometric impediment with 
I~ = (1/t c )T~ where the operator X~ is the integral op- 
erator on cf>, t c is the characteristic time of scattering on 
cuts or half-disks which depends on the concentration of 
scatterers n c . The power s = —1/2 corresponds to scat- 
tering on impurities/scatterers with fixed density in space. 
We note that s = 3/2 may correspond to a case of charged 
non-screened impurities distributed in the system plane. 
Eq. (]28[) can be transformed to the form: 



e 3 m 



2tt 

' {1+L0 2 T 2 ) 



°T c e(l+UJ 2 T 2 ) 



{(2 



- 2s)a 



(29) 



where v = vu, aijk =< UiUjT~Uk > and u is a uni- 
tary vector. The stationary current ji = aijkE u jE* k un- 
der linear polarized field is given by the following compo- 
nents of photogalvanic tensor (see Figs. I1I2[) : a xxx , a xyy 



and Re(a yxy ). At T 
expressions for these components read 

e 3 V F T 3 



(a degenerate Fermi gas) the 

(30) 



(2 - 2s) + 



7TT C (1 + U) 2 T 2 ) 

S {l-UJ 2 T 2 ) 



(1 + 

e 3 V F r 3 



(2 ^S^)d X xx 

(2 



7TT C (1 + W 2 T 2 ) 

S {l-U 2 T 2 ) 



(1 + U ¥) uxyv 
e 3 V F r 3 



2s)a x yy + 



7TT C (1 + U1 2 T 2 ) 

s(l-w 2 r 2 ) 



(31) 



(32) 



(a x 



yy 



Here, r, r c are relaxation times taken at e = ep. Now it is 
necessary to calculate two quantities a xxx and a xyy . They 
depend on the model of asymmetric scatterers. For the 
case of cuts the scattering probability has the form: 



W(<f/,tf>) 
1 



TVT^COS^ 9(cOS(f>')S((f)' + 4>-w) 



cj>' cos <fi #(cos <p)0(— cos <f>')] . 

Using Eq. (l33|) we obtain a xxx — 1/48, a xyy 
a result we have in this model: 



e 3 V F r 3 



48ttt c (1 + u 2 t 2 ) 2 
x[2-s+(2-3sV 2 r 2 ] , 
e 3 V F T 3 



x[2 



' 48ttt c (1 + lo 2 t 2 ) 2 
s + (2 - hs)uj 2 T 2 ] , 
e 3 V F r 3 



24ttt c (1 + uj 2 t 2 ) 2 
x[(3-2s) + (3-4s)w 2 t 2 )] . 



(33) 
-1/16. As 

(34) 
(35) 
(36) 



The formulas (|34I35I36|) with s = —1/2 follow also from 
the exact solution of the problem for cuts obtained in . 

According to the previous Section for the semidisks 
model the scattering probability is 



w(<p',4>) 
i, ,■ 



cos <p' 0(cos <t>')8{4>' + 4> -it) 

^)\[e{j>-ci>')e(-<j>'-<t>) 

(37) 



This probability W leads to the following results for the 
photogalvanic coefficients in the semidisks model: 



-a x 
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Fig. 5. Dependence of the rescaled velocity of ratchet Vf /Vf 
on the rescaled collision time r c /ri obtained by numerical sim- 
ulations of the cuts mode (symbols) for u)Ti = 4.7; 7.05; 9.4 
(red, black, green symbols from top to bottom respectively). 
The full curves show the theoretical dependence Eq. (|43l45p for 
corresponding uin multiplied by a numerical factor Q = 0.54. 
Here, the equilibrium Fermi-Dirac distribution has the tem- 
perature T/Ep = 0.1 and eEriVp/Ep = 4.2; the polarization 
angle 8 = 0. 



e 3 V F T 3 



127TT C (1 + 0J 2 T 2 ) 2 



[2 - s + (2 - 3s)lu 2 t 2 } . (38) 



In the numerical simulations for the cuts model it was 
taken that t is independent of energy (s = 0). In such a 
case from the above Eqs. (|34l39B6p and Eq. (13"8"1> we obtain 



1 

— tie{a yX y) 



\/2r, 



247n A /mepr c (l + lo 2 t 2 ) 

&xxx — &xyy — Rc^y^y) 



(cuts); 



(39) 



*xyy 
e 3 T 



3 - 3 V2n e 



6r c m^/m£i?(l + uj 2 t 2 ) 



(semidisks), (40) 



where n e is the electron density related to the current by 
the relation j = en e Vf . 

Let us also note that at r = const for the case of 
Maxwell equilibrium we obtain 



a xyy — ^R- e ( a yxy) 
3 3 



2-nn e T i e 



24mVmTT c (l + lu 2 t 2 ) 

x ^xyy f^K^yx?/) 

\/2nn e T 3 e 3 
6mV toTt c (1 + lo 2 t 2 ) 



(cuts) , 



(41) 



(semidisks) . (42) 



It is important to note that the relation a xxx = —a xyy , 
which is valid for semidisk antidots at any s, implies that 



the total photogalvanic current is zero for depolarized ra- 
diation. For a general form of scatterers and r dependent 
on energy (s ^ 0) even depolarized radiation produce 
nonzero photogalvanic current. This is for example the 
case for the cuts model with fixed density of impurities 
where s = —1/2. 

We should again emphasize the difference between for- 
mulas obtained using the kinetic equation approach and 
simplified way which leads to Eqs.([6][7]). The difference 
in numerical factors appears due to different methods of 
averaging: in Eqs.([5][7|) it is assumed that the scattering 
events are randomly distributed in time, while in the ki- 
netic approach it is assumed that the scattering events 
are randomly distributed in space. A simple way to ob- 
is to consider the mean force applied 
gas to an asymmetric scatterer. This 



tain Eqs. (|39M2|) 
from the electron 
force is determined by the momentum production on the 
scatterer, or, in other words, residue of flows of momenta 
of incident J v x pD f (p)dp/2TT 2 and scattered particles. If 
the anisotropy is weak, the momentum production can be 
found substituting the expression for distribution func- 
tion (JTHJ) into the collision integral The flow velocity 
is determined by equating this force per unit area to the 
friction of electron system, — P/r, where P = mwfji e is 
the full mean momentum of electron gas. The result ob- 
tained in this way is the same as those obtained from the 
kinetic approach if the mean free time does not depend on 
the energy (Eq. [391132]) . 

For a direct comparison with numerical simulations 
it is convenient to rewrite Eqs. (|39I40[) to obtain explicit 
expressions for the ratchet velocity. For the Fermi-Dirac 
distribution at T <C Ep this gives for the cuts model 



Vf = — 

24 



- cos(2<9) \ 
3sin(20) J 



and for the semidisks model 



where 



Vf 



5' 



_ Sf- 008(20) \ 
3in(20) J 



(43) 



(44) 



e 2 T 3 V2E 2 



m\J mEpT c (\ + uj 2 t 2 ) 



(eETV F ) 2 T 
2E 2 f t c (1 + U 2 T 2 



(45) 



To compare the polarization dependence given by 
Eqs. (|43I44I45P (see also Eqs. (|41I42|) ) with the numeri- 
cal results shown in Figs. 13141 we should take into account 
that in the numerical simulations of the cuts model the 
collisions with cuts take place always after time t c (ran- 
dom positions of cuts in time) while the computations 
above assume fixed distribution of scatterers in space that 
leads to a factor 3 sin(20) in Eq. ((431 instead of 2 sin(20) 
in Eqs. (|9ll9p . A part of this, the functional dependence 
on parameters is the same in both computations. A slight 
difference in a numerical coefficient related to different 
ways of averaging for cuts randomly distributed in time 
(Eq. HH)) and in space (Eq. (03}). 
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The comparison of the theoretical dependence on r, r c 
(Eqs. (|43I45|) with the numerical data in the cuts model is 
shown in Fig. [5] at different values of microwave frequency 
and 8 = 0. It shows that the numerical data are well de- 
scribed by the theoretical relation Vf = —QS/24 with a 
numerical factor Q ss 0.5. This confirms that the theory 
gives a good description of the numerical results. 

It is also interesting to compare the theory for the disks 
model with the numerical data of [28] (see Fig. 8 there). 
Using the value of t c from Eq. ([1]) for a low density at 
R/r d fti 4 and eEr d /E F = 0.5 we obtain from Eqs. (1441451) 
Vf /Vf ~ 0.14 instead of the numerical value Vf /Vp ~ 0.05 
(according to Fig. 8 in [28] at 8 — and R/r d = 4we have 
lot 0.8 since Tj S> t c ). This means that a numerical 
factor Q between the theory Eqs. (|44l45p and the numer- 
ical simulations [28] is Q » 1/3. Such a correction can be 
considered as rather satisfactory since the theory assumes 
low density approximation which is not yet well justified 
at R/rd = 4. Also the theory of kinetic equation works in 
one-collision approximation which may work not so well 
for completely dynamical systems of semidisks lattice. 

Finally, let us return to a simple derivation of the pho- 
togalvanic current given in a paragraph after Eq. ([42]) . Ac- 
cording to this consideration the directed flow appears on 
a time scale of the order of r. For a purely dynamical sys- 
tem, like the semidisks Galton board, we have Tj = oo and 
t = t c . We note that r c is given by Eq. ([TJ and is propor- 
tional to the inverse Lyapunov exponent of chaotic dynam- 
ics in absence of microwave driving. This leads to an inter- 
esting fundamental question about the photogalvanic cur- 
rent in a purely dynamical system. Indeed, in absence of 
friction the semidisks Galton board with a monochromatic 
driving is a Hamiltonian system which has time reversibil- 
ity property (like cos cut), while the appearance of directed 
current on a time scale r c breaks time reversibility. The 
dependence of the photogalvanic current on parameters of 
a real relaxation process is clear from the phenomenolog- 
ical point of view: the current changes its sign with the 
time inversion while the electric field squared keeps its 
sign (Eq. (|21jl ). This means that the photogalvanic coeffi- 
cient a in Eq. (|2ip should contain the relaxation constant. 
By its own, the kinetic equation is irreversible even if only 
the static potential scatterers are taken into account. In 
fact, a purely dynamical system is reversible and hence the 
stationary current should vanish for it. In principle this 
point can be explained by the dynamical chaos where the 
time reversibility is broken in practice for a coarse-grained 
distribution. However, a more delicate point is the ques- 
tion about the detailed balance principle. In a dynamical 
system it means that the transition probabilities are pro- 
portional to a measure in the Hamiltonian phase space 
and if all phase space is chaotic there should be no global 
directed current on large time scales. In a sense the nu- 
merical results |26j are in favor of this statement since 
there a small friction force F f = — 7P gives the velocity 
of stationary flow Vf ~ jR which disappears in the limit 
of 7 — > 0. However, in the limit of small 7 the situation 
is somewhat specific since the average steady state energy 
E s , analogous to temperature, grows with a decrease of 7 



as E s ~ ((eE) 2 /j) 2 / 3 ~ T that leads finally to the rela- 
tions similar to those given by Eqs. (|9|15|41|42|) . The nu- 
merical simulations performed in |27j have been done with 
the Nose-Hoover dynamics (see [30]) which can be viewed 
as a purely dynamical time reversible Hamiltonian sys- 
tem in an extended phase space with additional variables. 
Thus the photogalvanic current can appear in a Hamil- 
tonian asymmetric system with monochromatic driving. 
However, the directed flow we discussed takes place only 
in some part of total phase space corresponding to phys- 
ical variables and it is possible that the total current in 
total extended phase space still remains zero. Indeed, in 
principle it is known that in chaotic Hamiltonian systems 
there may be two separate components (e.g. one with a 
regular motion and another with a chaotic motion) with 
a directed current in each component but with the total 
current equal to zero (see examples in [31]). The ratchet 
analyzed in [27] with the Nose-Hoover dynamics can be 
such a case. Also, the contradiction can be resolved if to 
assume an existence of a certain time scale after which the 
photogalvanic current is stabilized. Thus, from a practi- 
cal view point we may say that in a dynamical chaotic 
Hamiltonian system signatures of directed flow appear af- 
ter the time scale t c but since the current velocity Vf it 
proportional to a second power of weak field Vf cx {eE) 2 
a relatively long time t r cx l/(eE) 2 is needed to observe 
this current in presence of chaotic fluctuations and it is 
necessary that the steady state distribution in energy is 
established on a time scale which is shorter than t r . Fur- 
ther studies are required to understand more deeply the 
problem of time reversibility in the context of the photo- 
galvanic effect. 

4 Effects of a magnetic field for the semidisks 
Galton board 

For experiments on photogalvanic current in asymmetric 
nanostructures it is important to know what are the effects 
of a magnetic field B perpendicular to the 2DES plane on 
the strength of current and its directionality. An analytic 
solution of the kinetic equation becomes much more com- 
plicated compared to the cases considered above. This is 
especially the case when the Larmor radius Rl of electron 
motion becomes comparable with the size of asymmetric 
antidots. Therefore, the numerical simulations in this case 
become especially important. For the semidisks Galton 
board the effects of magnetic field have been studied in 
[28]. They clearly show that the ratchet current becomes 
quite weak when the Larmor radius Rl becomes smaller 
then the semidisk radius r d . This is rather natural from a 
physical view point since in this regime the scattering on 
semidiscs is suppressed. 

However, in the regime with r d /RL ~ 1 a relatively 
weak magnetic field can significantly affect the directional- 
ity of photogalvanic current. This is illustrated in Figs. 16171 
obtained by numerical simulations with the method de- 
scribed in [28] . The results of these Figures clearly show 
that a moderate magnetic field can change the direction 
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0.008 



-0.004 - 



-0.008 



Fig. 6. Dependence of rescaled ratchet velocity v x /Vf in 
a;— direction on the ratio ra/ Rl proportional to the magnetic 
field B perpendicular to the semidisks Galton board. The sys- 
tem parameters are R/r c i = 2.5, T/Ef = 0.1, eErd/EF = 0.3, 
lui/fd = 2-kVf / (urd) = 9. The bottom (top) red (blue) curve 
corresponds to the polarization angle 9 — (9 = 7r/2). The 
rescaled mean free path due to impurity scattering is l/ra = 
VFTi/rd = 45; the rescaled computation time is cot = 3.1 x 10 6 . 




Fig. 7. Dependence of ratchet current angle tp on a rescaled 
magnetic field r d /Rh oc B for the case of Fig.[6]at 6 = 0; curve 
with circles shows numerical data, angle is given in radians. 



of current almost on 180 degrees (Fig. [7]). The angular de- 
pendence of Fig. [7] is not sensitive to the microwave field 
strength and therefore is not related to the Lorentz force. 
We attribute the origin of this strong angular dependence 
to a significant change of scattering process in the regime 
when ra/Ri ~ 1 related to multiple collisions of electron 
with a semidisk. 

To make a more close link to possible experimental 
studies we note that for the electron density n e = 2.5 x 
10 11 cto -2 , an effective electron mass m = 0.067m e and the 
semidisk radius ra = OAfim we have for the parameters 




-0.004 - 



-0.008 



Fig. 8. Dependence of the ratchet rescaled velocity v x /Vf on 
the rescaled ratio rd/Rh oc B for parameters of Fig. [6] at 9 — 
and various values of rescaled microwave frequency uird/VF = 
0.335 (brown), 0.447 (orange), 0.648 (magenta), 0.693 (red), 
0.805 (cyzen), 0.939 (yellow), 1.12 (blue), 1.39 (violet) (curves 
from top to bottom at rd/Rh = 0.8). 



of Fig. [5] the following physical values: Ep 100K, Vp = 
2.2 x 10 7 cm/s, w/2tt = 60GHz, l u = 2ttV f /uj = 3.6/im 
and I — VpTi — 18/iTO (such value of / corresponds to 
mobility of about 2 x 10 6 c77i 2 /V^5), a magnetic field B = 
0.075T corresponds to the Larmor frequency u) L w SAGHz 
and rd/Rh = 0.4. For these parameters and data of Fig. [8] 
the microwave frequency changes from 29GHz to \2QGHz 
when the rescaled ratio changes from wfrf/Vp = 0.335 to 
1.39 {cur d /V F = 0.693 for Fig.©. 

In absence of magnetic field there is no dependence of 
current direction on the microwave frequency. On the con- 
trary, in the regime r^/Rh ~ 1 the directionality of flow 
can be also changed by changing ui as it is shown in Fig. [5] 
We also note that the velocity of flow in x direction re- 
mains the same with a change B — > —B due to symmetry 
reasons and we present data only for B > 0. 

This angular dependence becomes weaker when the 
mean free path I decreases due to decrease of impurity 
scattering time Tj (see Fig. [S]). The results of Fig. [5] also 
show that the ratchet effect disappears when I becomes 
smaller than size of asymmetric antidot. Of course, the 
current is absent when the semidisks are replaced by disks 
(in this case the numerical data, shown by dashed curve 
in Fig. [9l are on a level of statistical fluctuations). 

The angular dependence of the photogalvanic current 
B and u> found in this Section is rather nontrivial and fur- 
ther studies are required to obtain a detailed explanation 
for it. 



5 Discussion 

In this paper we developed a theory which determines 
the strength and directionality of the phogalvanic cur- 
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0.008 




-0.004 



-0.008 



Fig. 9. Dependence of the ratchet rescaled velocity v x /Vf on 
the rescaled ratio rd/ Rl oc B for parameters of Fig. [6] at 6 — 
and various values of rescaled mean free path l/rd = VFTi/ra ~ 
45, 22.5, 11.1, 4.5, 2.25, 1.1 (full curves from top to bottom at 
rd/Rh = 0.8; for the physical parameters given in the text I 
changes from 18/^m to 0.45/^m); dashed curve corresponds to 
the case when semidisks are replaced by disks and Ijrd = 45. 



deed, the theory developed above considers only the case 
of noninteracting particles. In reality electron-electron in- 
teraction may play an important role, as well as plasmon 
modes can also affect the photogalvanic effect. Also at 
small scales with R ~ lOOnm the level spacing between 
quantum levels inside one lattice cell becomes of the or- 
der of 100GHz (at n e ~ 3 x 10 11 ) and quantum effects 
may play an important role. Also at R < lOOnm the col- 
lision frequency enters in teraherz range Vf/R > ITHz. 
Thus such asymmetric nanostructures can be used as room 
temperature detectors of radiation in teraherz range. In 
analogy with the GORE-TEX material used in common 
life for one-way transport of water and air [32] such ar- 
tificial asymmetric nanostructures may be considered as 
"NANO-GORE-TEX" material. The further studies of the 
NANO-GORE-TEX properties may bring new interesting 
applications including high frequency detectors and sen- 
sors. 

We thank Z.D.Kvon for useful discussions. This work 
was supported in part by the ANR PNANO project MI- 
CONANO and (for MVE and LIM) by the grant of RFBR 
No 04-02-16398, Program for support of scientific schools 
of the Russian Federation No. 4500.2006 and INTAS No. 
03-51-6453. 



rent in artificial asymmetric nanostructures. The theoret- 
ical results are in good agreement with detailed numerical 
simulations performed in this work and in [28]. We also 
find from numerical simulations that the directionality of 
photogalvanic current is very sensitive to weak magnetic 
fields. 

A microwave field of E = IV/cm generates in a lat- 
tice of semidisks with R ~ 1/im a current of about 0.2 
nA per lattice row at electron density n e ss 2.5 x 10 11 . 
For a structure of 100/im we have currents of about 20 
nA. The photogalvanic effect has classical grounds and 
exists at room temperature if the mean free path remains 
larger (or comparable) than the size of asymmetric an- 
tidot. This is the case of the Lund experiment [25 a with 
I ~ R ~ lOOnm (the polarization there corresponds to 
8 = 7r/2, triangles are used instead of disks). For a fixed 
ratio R/rd a decrease of R ~ I by a factor 10 gives a drop 
of current by a factor 100, however, the induced voltage 
drops only by a factor 10 since a resistance also drops 10 
times. While the experiments I23|25j demonstrate the ex- 
istence of the effect there are still no experimental data to 
be compared with our theoretical and numerical results: 
in the experiments the strength of current is known but a 
microwave field strength acting on electrons is not well de- 
fined. Also till now only two polarization cases have been 
analyzed experimentally (9 — in [23 and 8 = tt/2 in 
[25 ) . The directionality of current in these experiments is 
in agreement with our results but a detailed experimental 
investigation of the polarization dependence is still highly 
desirable. 

All these considerations clearly show that further ex- 
perimental studies of photogalvanic currents in asymmet- 
ric nanostructures are very interesting and important. In- 
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